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Abstract 

We consider the AdS4 RN black hole and work out the momentum dependent hydrodynamic 
analysis for the vector modes. We also perform the spectral function calculation of the dual field 
theory. As an application, we consider the permittivity and permeability and find that for low 
frequency regime, the index of refraction is found to be negative, supporting the claim made in 
ref |[20l for AdSs. We also find that at static limit the medium has the zero permeability, a character 
of the superconductivity . 
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The gauge/gravity duality [[T]-l3l has been proved to be the most fruitful idea in recent string theory 
and has been widely used in several different systems: from the quark-gluon plasma Il6]-fr3l and to 
condensed matter systems lfT4l[T6l[T7l . The super conductivity ifTSl and transport properties ifTSl were 
calculated for the 2+1 dimensional system near the critical point and dual background of the non-trivial 
scaling properties were also suggested lfT9l . In a very recent paper [[201, a new connection between AdS 
gravity and optics we suggested: the negative refractive index of some meta-materials can be described 
by the holographic duality. 

The most basic configuration is the asymptotically anti-de Sitter black hole with electric charge, 
which describe the finite temperature and chemical potential of the dual field theory. For the trans- 
port coefficient, one only needs to calculate the Green functions for the zero momentum case and 
corresponding analysis were done in lfT6l[T7l . On the other hand, the permeability request momentum 
dependent hydrodynamic analysis even for the small frequency limit. 

In this paper, we first perform the momentum dependent hydrodynamic analysis for four-dimensional 
Reissner-Nordstrom-Anti-de-Sitter black holes. The decoupled equations of motion and the Green 
functions have been worked out explicitly. We then study the spectral functions for the vector modes 
and as an application, refractive index for a strongly coupled medium were calculated. 
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In lfT6l[T7l . hydrodynamic analysis for the dyonic black hole for zero momentum case but with 
magnetic charge was performed. The zero momentum limit of our result is different from their result 
by taking the zero magnetic limit. This is because the limit of zero magnetic charge and that of zero 
momentum is not commuting. It is singular as was argued by the authors of ref. [[T5l . However our 
result gives consistent result for the conductivity in the literature ll22l[23l . We also calculated the 
spectral functions, permittivity and permeability and for low frequency regime, the index of refraction 
is found to be negative supporting the claim made in ref. [|20ll for AdSs. After finishing this work, 
we are informed that propagator structure with general quasi-normal modes has been studied in a very 
recent paper [|28l and hydrodynamic analysis at zero temperature is described in [|29ll . 



2 Basic setup 

The action for Einstein-Maxwell theory with a negative cosmological constant A = — is give by 
S=^ [ d^x^{R -2A- L'F^.F^') + \ f (i'x^Q + 4f / d'^^f^^ (2-1) 



J ^4 J ^4"^ 

where F^^, is the field strength of Maxwell field A^,, the second integral contains the Gibbons -Hawking 
term and the third one contains the counter term for the regularization. Here 7^1, is the induced metric 
on the boundary and 6 = 'y^'^Q/j,^ is the trace of the extrinsic curvature = — (V^^ni. + V^n^) with 
n the unit normal vector to the boundary directed outward. 

The 4-dimensional Reissner- Nordstrom- Anti de Sitter black hole is the solution of Einstein equation 
and the metric is 

ds' = '-^{-f{z)dt' + dx' + dy') + L^ (2.2) 
^ j[z) 

= -Qaz. (2.3) 

The function f{z) has the form 

/(z) = 1 + QV-(1 + Q2)z3. (2.4) 
The Hawking temperature of the black hole reads 

r^P-/)°. (2.5) 

Air 

We consider the linearized perturbations htx, hxy and of the Einstein equation and the Maxwell 
equation, g^,^ = gf} + h^^y, Afj, = A^^^ + A^. From the Maxwell equation we find 

fifA'J + u'Ax - QfhT - klfAx = 0. (2.6) 
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The linearized equations from the Einstein equation yield 

-hf + -hf + ^Qz^K + jK + = 0, (2.7) 



K" + ^-zH-K + - {k,h^, + ujhl) = 0. (2.9) 



w/if - 4wg2M^ + kyfh''y' = 0, (2.8) 

In these equations we have transformed the variables — )■ A^a, cu — )■ coa and A; — j- A;a, so that 
all the coordinate, momentum as well as field variables are dimensionless. Therefore that and hy 
were defined to be dimensionless from the beginning. At the end of the calculation, we can restore the 
original variables by reversing the scaling. 

In order to have a decoupled equation, we define a new variable 

^± = %-AQA, + ^A,. (2.10) 
Then equations (|2.6I) . (|2.7I) and (12.81) reduce to two decoupled equations for $± 



*± + H^*± + T^-^i - ^F^-^i = 0, (2.11) 



where C± = V9(q^+i)^+i6.^ q^_ 

3 Hydrodynamic analysis 

We now seek a transformation 

<!>{z) = giz)E{z) (3.1) 

such that the coefficient of ^(z) goes to zero when we take the hydrodynamic limit, uj,k — )■ 0. Here $ 
represent one of $±. 

k^ + C^Qz , 2g'{z) , f{z)g'{z) , ^"(^)^ , ^ , , .g'i^) ^ , 



f{zy fiz) zg{z) g{z)f{z) g{z) J ' ' \z f{z) g{z) 

(3.2) 

We need to find g such that S's coefficient in eq. (l3.2l) is of order k"^ or higher in hydrodynamic limit. 
The constant C± can be expanded in a series of k'^ 

C^^'JllJll^^^^Om (3.3) 



Q 3(l + g 
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3.1 For 

We request that g{z) satisfies following equation 



^^^^^^)%(.)+r^ + ^l.'(.)./(.)^0. (3.5) 



9i^) = - - (3-6) 



The solution for (|3.5I) regular at 2; = 1 is given by 

1 4g2 

^ 3(1 + Q 

where we used the linearity of equation of motion to normalize the solution. 

By imposing the infalling boundary condition and taking out the near horizon limit, 

S(z) = (l-z)^F(z), (3.7) 

withi/ = -icj/47rT.Q 



+f-(,)(_2^ + m + ^M£))+r'W^0 (3.9) 

\\- z f{z) z g{z) J 

In the long- wavelength, low-frequency limit, we can expand F(2;) in a double series with respect to to 
and k 

F{z) = Co + coFi{z) + eF2{z) + 0{u^e) (3.10) 



, 1 f ^Co 2iCo zCof'jz) 2iC,g'{z) \ 

^^-Q-\{-l + zf^ z-z^^ f{z){l-z)^ g{z){l-z)) 

^ [ ( (3-g4-i) + wMw) ^'^"''^^^ 

^ _ Co r jl x\?, + Q\?> + Ax)) g{xf dx 



3(1 + Q2) 



j-x- (3 + g-(3 + 4x))^(x)-cfx ^ 
/ ^ s . ^2 dy:=CoJ{z) (3.14) 

Ji y f{y)9{y) 



'Here both w and T are scaled by a to be dimensionless: T = (3 — Q^)/47r and = —iuj/{3 — Q^). 
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Later we will need following results for H{z) and J{z): 



H{z) = H{0) + ifi^ + ■ ■ ■ , J{z) = J(0) + Jiz + --- 

with Hi = 7T, (3.15) 

9(3-g2)(i + g2)2' 

27 + 63g2 + 29g4 + gg^ ^^^^^ 

Ji = n . (3. Id) 

27(i + g2)=^ 



3.2 For $ 



Since C_ ~ 0{k'^), we do not need transformation factor and we set g{z) = 1, unlike $+ case. By 
doing as before, we get 

u{-2 + z + zu) uf'{z) (^'-^'/W (l-3(W))) 



F'{z) + F"{z) = (3.17) 



z{z-iy {z-l)f{z) f{z^^ 

'2 , f'{z) , 2z/ 



,z /(z) z-1 

Expanding F{z) = Co + ujFi{z) + k'^F2{z) + 0{u^, k"^), we have 

^ ^ +'- + m) + (-^Q^]F.'iz)+FAz) ^ 0, (3.18) 



(3-g2)(i-z) vi-^ ^ fiz)j \Z fiz) 

)+i^2(^)(y|| + ;)+i^2'W = (3.19) 



Co AQ'z 



f{z)\ 3(l + g2 
and their solutions are 



F2 = -TTTT^fl--'! (3.21) 



3(i + g2)V z 

respectively. 

Let us fix the constants Co and Cq by imposing boundary conditions as 

lim/i^(z) = h^, \imhl{z) = h^, limA^ = (3.22) 

z— >-0 z— >-0 z— >0 

It would be simple by taking a derivative of $±(-2) and using equation (12.71) . This gives the relation 

^2$; - C±zA'^ = (^Ph^ + ukhl^ - C±A^. (3.23) 
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In order to determine Cq and Co, we need to examine the above equation at the boundary 

lim(22$'^ - C±zA'^) = (^'^^ + ukff^ - C±A^ := L± 

We have a few remarks on eq. (13.241) . 

• Only i singularity of $ near boundary can contribute to the left hand side. 

• The boundary value of zA^ is 0. 

• The equation ( 13.241) is correct even for the finite k and tu. 

Let 

$± = — + ^± + T±log^••• 
be the expansion near the boundary. Then 

<!>+ = Co (1 + iujH{0) + A;V(0)) = -L+, 



and 



n_ = Co + o{kuj,u^) 

T± = 

where Hi, Ji are given in eq. (|3.16|) and eq. (13.161) respectively. 

The constants Co and Cq can be read off from the <I-|- and L± parts, respectively 



Cn 



Cn 



1 + iuHiQ) + F J(0) 
ehf + ujkhi^+^^^A, 



3(1+Q2) 

where -^/(-z) and J{z) were defined at eq. (13.121) and eq. (13.141) respectively. 
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Let us turn to the A'^{z) and hf near the boundary. From the definition of the master fields, we can 
solve hf and Ax{z) in terms of the $-1-: 



$1 - 



C+-C. 

For example A'(e) can be calculated from the observation: 

-n 



Afr, — Z 



C 4. - c_ 



(3.33) 



A'(e) 



C<-C. 



(3.34) 



where H-i- are defined in eq. (13.251) . Now, we can evaluate the solutions for the fields near the boundary 

z = e: 



hf 



hi' 



a' 



( 



a2 \ 



kuhl + u^frl ) - 



a iuj — Dk"^ 

e2 kulvl + u'^hl + AQDkojA^ 
a iu — Dk"^ 



A'{e) = -QDa 

+ A^i iu 



k'^h^ + ukh^ + AiuQA^ 



iuj - Dk^ 



AQD^ak^ {Klf 



9a(l + g2)2 9(1 + Q2)3^2 



k' + 0{e' 



(3.35) 
(3.36) 

(3.37) 
(3.38) 



Notice that a factor were restored and D is the diffusion constant 

1 



D 



3a(l + g2 



(3.39) 



For the Green functions, we need on-shell action that gives finite and quadratic function of the 
boundary values. Using the fact 



5S= d, 



dL 



and deleting all the contact terms, it is given by 



dz + / E.0.M.5(l)'dz, 



(3.40) 



r d^k 



^ ( h^hr - f{z)hlhl' ] - f{z)AA + QaAh- 



+ 



(3.41) 



2-^0 



The factor | comes from —3/4+1 where —3/4 is from the original action while +1 is from the Gibbons 
Hawking term. 
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Following ref. [4J, we discard the contribution of the horizon as a prescription for computing ther- 
mal Green's functions in Minkowski space. From the definition, Gi-j = —5'^S/5(pi5(pj, one finds the 
correlators in the hydrodynamic approximation 



_ ^2^2 ^2 

Gxt;xt 7i 9" ' ~ 7-1,9) Gxyxy ~Z 9" ' ~. rTTo' (3-42) 

2fi:| ^a; - Dk^ 



G,,„ = i!^ f-^!^ - 4QZ)V) , (3.44) 

_ 2^ /4za^aQ2D _ . - Q^f (9 + 2ig^-Q^+3g6)\ 

The transport coefficients can be read off from the green function. The DC conductivity a is given 

by 

a = - lim -ImG,.,.(a;, A; = 0) = (3.47) 

Notice that this is dimensionless and vanishes as square of the temperature ( ~ T^) near the zero 
temperature while it is constant in high temperature. 



4 Spectral function 



To calculate spectral function of current-current operator, we need to know the boundary value of the 
fields and theirs conjugate momentum. For the numerical computation, it is better to introduce new 
master variable ^f-t in this section which is different from ^± by a factor 1/z, 

^± = = -K + (C± - AQz) Ax. (4.1) 
z 

We remark that the new master field has nothing to do with \1/ in eq.(3.1). The equation of motion 
in terms of the new master variables are 

7** ^ 7 [~f ^^'^^ ~ ~ 

where C± is 



+ T*± + - I - <^±<52 - - ) *± = 0, (4.2) 



^ 3 + 3g^ ± v/9(Q2 + 1)2 + 3(i + g^/ A , ( 4gfc , , ,,,, 
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From the equations for master field we can get the spectral function of master fields. We however 
need the spectral function of original variables not the master fields itself. So let us first find the series 
solution of /if, Ax which defines the conjugate momentums nh, it a 



ht = h^t - -Ziz"^ + -ifhZ^ + 



2 1 
2 3 

Ax = Ax + TiaZ + -- - , (4.4) 



where Zi = kh^ + tuhy. The coefficient (3i is computed from eq. (12.71) : The master variables have 
series solution near the boundary, 

^± = ^± + n±z + ■ ■ ■ . (4.5) 
The eq. (14.11) can be written as a matrix form [[2111 

Similarly, if we introduce the transformation matrix Rq , 

Ro = I I 1 , (4.7) 

the boundary values and of master fields and those of original fields are related by Rq. 

we use the linear response relation II-i- = ^±G± in eq. (14.81) . Then the conjugate momentum tt/j, ifa are 
written as 

r ) " ( ) ' ^^^'^ ^ ^ Ro 'Diag(M/+, 

^ ^ 1 ( c4kz,-c^Ax) c4-kz, + c.Ax)\ 

y -kZi + C+Ax kZi - C^Ax ) 

From the boundary action eq. (l3.41l) . we compute two point function in terms of boundary values 
(/if, Ax) and conjugate momentum of master field n±: 

6h^6hf 6h^dh^ SAx6Ax 
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ImG, 




Figure 1: Spectral function of xt,xt component, Im Qxtxt- Left : with fixed k=0.5 varying Q=0.5(thin), 
l(thick), 1.5(dashed). Right : with fixed Q=l, varying k=0.2(thin), 0.5(thick), l(dashed) 



The two point function for and is related by Ward identity, therefore the correlation functions for 
each components are 



^2q,3 t'xtjXt ^ ^ ; j2 ^r^^ xy ,xy ^xy,xy ^2 ^^'*>2^i 



i^xt,x yx,xt 1^ 



70 ^ XL,X Z/X,XL ^ 

ra-^ C+ — C 

h 



<^xx = yx,x = — 7; • (4.11) 



4.1 holographic recipe for Green function 

In general, the second order differential equation has two independent solution and these two indepen- 
dent solution have two different integration constant. We represent the near the boundary solution by 
$j and near horizon solution by 0j. Then, 

$1 = + </>i = (l-n)--(l + -.-) 

$2 = + ■■■), 02 = (l-w)"(l + ---) (4.12) 

where u = is the UV boundary or AdS boundary and u=l is IR boundary or black hole horizon when 
we consider finite temperature system. Near the boundary, u=0, $i and $2 are the local solutions 
and A± is the solution of indicial equation near the boundary, A+ > A_ where A+ is the conformal 
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Figure 2: Spectral function of xx component, ImQ^x- Left : with fixed k=0.5 varying Q=0.5(thin), 
l(thick), 1.5(dashed). Right : with fixed Q=l, varying k=0.2(thin), 0.5(thick), l(dashed) 



dimension of an operator and A_ is the dimension of the dual source field. Near the horizon, u=l, 
(pi , (p2 are the local solution and they are identified with infalling and outgoing solution respectively. 



A{rx>, q)$i + B{rx>, q)<l>2 = C{vo, q)0i + V{rx>, q)02 
01 (m) =^(t»,q)$i(M) + i3(i»,q)$2(M). 



(4.13) 



where A{vo, q) and B{rx>, q) are two integration constants and will be identified with the boundary value 
of the bulk field and its conjugate variable, respectively. To get the Green function it is first to evaluate 
the on-shell action. 



'S'on-shell = ["^ j d'^x^g''^ EaEp'\u-. 



Um / drxA 



^-{A++A_-l)^2 



A_-l 



A+-1 



^-A 



(4.14) 



The first term in the last line is obviously divergent and it should be renormalized by introducing 
counter term. The second term is Green function which is obtained by differentiating on-shell action 
with boundary value. 

More specifically, the local Frobenius solutions of our master variable near the boundary is 



^±(w, A;)^^(z, k) + n±(a;, k)-^1{z, k) 



1 + aiz + + 



^l{z,u,k) = z{l + (3iz + (32Z^ + ■■■) (4.15) 
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Figure 3: Comparison between numerically (solid) and analytically computed in hydrodynamic limit 
(dashed) spectral function with Q=0.5. Left is for G°(w) and right for ^^(w). These two functions are 
well matched in small w region but deviated with each other. 



What we want to calculate is Green function of each field and it is identified with conjugate momentum 
of boundary field dH. In eq. (14.141) . A_=0 for our master variable, so that kind of divergence is absent. 
Then 

fl±{u,k) d^(j)i{w,k,z) 



(4.16) 

2=0 



i!±{u,k) (j)i{w,k,z) 

Eq. (14.111) and eq. (|4.16l) give us Green functions of original fields. Only remaining is to compute II-i- 
numerically. 

4.2 small k expansion 

In this subsection, we will show how to compute G^xl with our previous master field. From eq. (14.111) . 
the expression for current-current Green function is 

c+g+ - C-G- 

yxx = — c - C (4-1 7) 

and Gxx is expanded near k=0 as G^x = C^^^ + k'^G^'^^ + ■ ■ • , so 

g = g(o) + k^gi^) + . . . (4.18) 

where = g+{w, k = 0) because (7_(0) = and ^(2) is 

G^'^ = I ^+(0) - ^40) 1 + G'JO). (4.19) 



c+(o) 
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To find second order solution for Q±, we also need to expand Q± itself as 

+ ■ ■ ■ = g^j;\w) + e g^^\w) (4.20) 



y±[w, k) = y±{w, k = 0) + k 



dk^ 



k=0 



First we find the infalling solution of the master variable which we call (p^ ' and expand it in momentum, 
k: 

= + ^ . . _ ^(0) ^ ^2^(2) ^ . . . _ (4_21) 

we then plug it into the eq. 14. II and expand the equation of motion. The zeroth order equation is 

^f" + ^j^f + j{j-Q' cm - j) = (4-22) 
and second order equation is 



+ ^^'^> + y (^y - C±(0) - ^1 = 0. (4.23) 



Numerically, , Q^^ are computed as 



ef = l.m^. 6f = nm^. (4.24) 



5 Negative index of refraction 

At the beginning, let us first review the electric permittivity e(uj) and the magnetic permeability h(uj) 
and their dependence on the Green functions. At the leading order of the electromagnetic coupling, the 
optical properties of the medium can be described by the linear response to an external electromagnetic 
field. In an isotropic medium with spatial dispersion, the electric permittivity e(uj, k) and the magnetic 
permeability /i(a;, k) depend on both the frequency and the wave vector k. The refractive index for the 
transverse modes has the form 

n'^{u,k) = e{LL!,k)fi{u,k) (5.1) 

Without the dissipation, e(uj, k) and fi{uj, k) are real. However, in case of that the medium is dissipative, 
e(a;, k) and /i(a;, k) will have imaginary part. So that the refractive index n is a complex quantity. The 
real part of n is the refraction index while the imaginary part encodes the information of dissipation. 
The refraction can be negative if Re{e) and Re{fi) are not simultaneously negative. For the dissipative 
case, the negative refractive index is equivalent to require 

riDL = \e{co)\Reifxiuj)) + |/i(u;)|i?e(e(a;)) < 0. (5.2) 
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In the Laudau-Lifshitz approach to electrodynamics of continuous media [|24] - I26l . the transverse part 
of the dielectric tensor is determined by the transverse retarded Green function as follows [|26ll27l 

er = l -GT{u,k), (5.3) 

where is the 3-dimensional electromagnetic coupling and Gt{uj, k) is the transverse part of the Green 
function. 

By expanding the Green function as Gt{uj, k) = G^p{uj) + k'^G^^\uj), 

(,^_ 9 + 2ig^-g^ + 3g6 AiQ^ 

9(i + g2)3 9(i + g2)2w; ^ 

the electric permittivity and effective magnetic permeability are given by 

e(w) = 1 —G^^> {co, k) = l + I ^ (5.5) 

id U 

= 3— . (5.6) 

l + ATie^GPicu) 

we use the eq. (I3.47|) for the expression of e(w). For sufficiently small value of e, the eq. (15.61) can be 
expanded as 

= 1 - Atic^gP (w) + C(e^) (5.7) 
From the above equations, we can see that at low enough frequencies, udl can be negative. From eq. 
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Im fi 




0.005 0.010 0.015 0.020 0.025 0.030 




0.005 0.010 0.015 0.020 0.025 0.030 



Figure 5: Real(left) and imaginary(right) part of transverse magnetic permeability fi{w) with e =0.5 and 
Q=0.5. 



(15.41) . (15.51) and (15.61) . in hydrodynamic limit, we find that 

AweH^a ft (3-g2)2 i 



1 + 



1 + 



AneH^a f9 + 21Q^ - + 3Q^ AiQ^ 



9(1 + Q2) 



9(1 + Q2)2y; 



w 9(1 + (y^)a 



(5.8) 
(5.9) 
(5.10) 



The last equality can not be valid for small enough frequency, causing the difference of our figure for the 
permeability from that in the ref. [|20ll . Our result suggest that at static limit, the medium exhibiting 
the character of the super conductivity, i.e, the zero permeability. 

It is clear that the value of n^L = |e(c<;)|i?e(/i(c<;)) + \fi{uj)\Re{e{uj)) can be negative at the low 
frequency regime, since the real part of e(a;) acquires negative value, while the imaginary part of /i(a;) 
is positive. 



6 Conclusion 

In this paper, we worked out the momentum dependent hydrodynamic analysis for the vector modes of 
charged AdS4 black hole. We also calculated the spectral function of the dual field theory and as an 
application, we calculated the permittivity and permeability and found that for low frequency regime, 
the index of refraction is found to be negative, which support the claim made in ref. [|20l for AdSs. 
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It will be very interesting if one can extend our analysis to the sound modes, which will enable us 
to calculate the longitudinal permittivity. It will be also interesting to workout the the case involving 
the magnetic charge. 
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